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MARKED-UP OF THE TRANSLATED INTERNATIONAL APPLICATION 

[Descriptionl TITLE OF THE INVENTION 

COMPUTER-AIDED METHOD FOR PARALLEL CALCULATION OF THE OPERATING POINT 

OF ELECTRICAL CIRCUITS 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to German Application No. 199 
27 301 .4 filed on June 15. 1999 in Germany, and PCT Application No. PCT/DEOO/01 7754 filed 
on May 30. 2000. the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The computer-aided simulation of electrical circuits has attained increasing importance 
in the deyelopment of yery large circuits, that is to say circuits haying a yery large number of 
elements. It is particularly in the deyelopment of computer chips haying a multiplicity, for 
example seyeral hundred thousand transistors, that serial processing for the determination of 
the circuit quantities by a computer has been found to be unusable because of the excessiye 
time consumption. 

[0003] In WO 98/24039. therefore, it is proposed to partition a large circuit and to haye the 
partitions processed by different computers in each case. 

[0004] In the calculation, the operating point, that is to say the potentials of all nodes is usually 
determined first as basis for further analyses such as. for example, transient or alternating- 
current analyses. 

[0005] For the parallel calculation, an implementation of the Newton method is proposed in 
U. Weyer, Q. Zheng et al.: "Domain Decomposition Methods for Circuit Simulation" 
(Proceedings of the 8th Workshop on Parallel and Distributed Simulation. PADS '94 Edinburgh. 
Scotland, UK, pp. 183-186, July 1994) and in U. Weyer. Q. Zheng: "Parallel Transient Analysis 
for Circuit Simulation" (Proceedings of the 29th Annual Hawaii International Conference on 
System Sciences, pp. 442-447. 1996). The disadyantage is that conyergence can only be 
achieyed here when sufficiently good estimates of the operating point are ayailable, due to poor 
conyergence characteristics. As a rule, howeyer, such good estimates can be achieyed with 
difficulty or not at all in the case of large circuits. 

SUMMARY OF THE INVENTION 
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[0006] [Thel One aspect of the present invention is, therefore, based on the object of creating a 
computer-aided method for parallel calculation of the operating point of electrical circuits which 
ensures simple, reliable and fast calculation of the operating point of the circuit. [According to 
the invention, this object is achieved by the method as claimed in claim 1 .] 

[0007] In the method, the charging method, which has very good convergence characteristics, 
is used for the parallel calculation of the operating point. 

[0008] As described, for example, in H. Spiro: "Simulation of integrated circuits" (2nd edition, R. 
Oldenbourg Verlag Munich, Vienna 1990), the charging method uses the dynamic elements 
such as capacitances and/or inductances existing in an electrical circuit for calculating the 
operating point via a pseudo-transient analysis. In this process, the dynamic elements are 
charged up step-by-step to a value which, in principle, can have any magnitude, particularly to 
the value "1", starting from a value zero. 

[0009] In an embodiment according to the invention, dynamic elements such as, for example, 
capacitances or inductances having predetermined known values can be additionally installed at 
suitable places, for example at or between a number of nodes. This makes it possible to avoid 
problems in the calculation due to the number of dynamic elements being too small. 

[0010] For this purpose, the dynamic elements additionally inserted into the circuit can be 
reduced step-by-step to zero from an initial value, for example tending to infinity, for which the 
operating point calculation is trivial, so that the original circuit is again simulated. 

[0011] In a preferred embodiment of the invention, at least one node, that is to say a junction of 
at least two current paths, of the circuit is connected by [means of] a capacitance to a 
predetermined value having a predetermined potential. 

[0012] In a further development of the invention, however, a capacitance, the second terminal 
of which is in each case connected to a predetermined potential, for example to ground, can 
also be connected to each node via all partitions. This procedure has the advantage that the 
calculation of the operating point, i.e. of the respective potentials of the individual nodes, is 
trivial for the circuit for an initial value for the capacitances which tend toward infinity, due to the 
equation, which is then explicit, being solved. Changing the value for the capacitances step-by- 
step then makes it possible to change the circuit simulation by suitable new selection of the 
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value for the capacitances until a calculation of the operating point of the circuit is obtained for a 
value of the capacitances tending toward zero or almost toward zero. 

[0013] By also implementing this charging method for parallel calculation of an electrical circuit, 
a very large circuit having a multiplicity of transistors can be advantageously calculated in a fast 
and simple manner even with a smaller number of dynamic elements. Various procedures are 
conceivable for the respective new determination of the value for the capacitance with the aim 
of allowing this to go toward zero, the decision criterion being the degree of difficulty of 
calculating the operating point of the respective preceding step. 

. ^ [0014] The operating point obtained by solving a nonlinear equation for a particular value of the 
Q capacitances in each case can be solved, for example, iteratively by [means of] the Newton 
'2 method. For choosing the next value for "C", the number of iteration steps necessary for the 
iQ solution can then be used until the value drops below a predetermined value for "C". [Further 
advantageous embodiments of the invention are obtained from the subclaims.] 

BRIEF DESCRIPTION OF THE DRAWINGS 

□ [0015] These and other objects and advantaoes of the present invention will become more 
lU apparent and more readilv appreciated from the following description of the preferred 
embodiments, taken in conjunction with the accompanvina drawings of which: 

[0016] In the text which follows, the invention will be explained by means of an exemplary 
embodiment shown in the drawing. In the drawing, the single figure shows a flowchart which 
represents the individual method steps of the method. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0017] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanving drawings, wherein like 
reference numerals refer to like elements throughout. 

[0018] According to the drawing, a simulation of an electrical circuit is the starting point 1 , for 
example in the circuit description language SPICE. This is partitioned in a first processing 
process 2, as described, for example, in WO 98/24039 as "Clustering method", so that individual 
partitions or parts of the circuits are obtained which can be calculated with the same degree of 
difficulty, if possible. 
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[0019] Naturally, other applicable variants of partitioning are conceivable such as, for example, 
the "ratio cut method" described in N. Frohlich, B. Riess. U. Wever, Q. Zheng: "A new approach 
for parallel simulation of VLSI circuits on a transistor level" (IEEE Transaction on Circuits and 
Systems - I: Fundamental Theory and Applications, Vol. 45, No. 6, June 1998, pages 601 to 
613), or even an arbitrary "manual" division. 

[0020] In a further process step 3, a grounded capacitance is added to each node of the circuit, 
that is to say to a junction of at least two conductors or current paths, respectively. Naturally, it is 
also conceivable to add a capacitance, the second terminal of which is connected to a 
predetermined potential, in which case both potentials and values for the capacitances of each 
node can differ. For reasons of efficiency, stability and economy of calculation, each capacitance 
is connected to the same potential with its second terminal, for example ground, an identically 
high value CO also being selected for all capacitances in step 4 for the same reasons. 

[0021] For this value CO, the operating point is then calculated for each partition or part-circuit 
in a further step 6, in which the required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being exchanged and included in the 
calculation of the operating points of adjacent partitions. 

[0022] In this process, a partition which will be called "master" in the text which follows can 
advantageously take over the control of the charging process for reasons of efficiency. The 
master then determines the initial value CO for the capacitances, and it is also conceivable to 
predetermine the initial value externally, for example by the user. This value CO is then 
transferred to all other partitions, called slaves in the text which follows. Following this, the 
operating point is calculated both in the master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the energy sources existing in the 
circuit, is present at the circuit. 

[0023] Starting from an overall problem to be solved 



where m is the unknowns of the master, 

p is the number of partitions, 
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Si is the unknowns of partition i, 
the charging process leads to the system of differential equations, 

f(x,t) + D^ = 0 



where 

D = diag (C, ... C, 0, 0) and 
t is the time. 

[0024] This system of differential equations can be solved, for example, by [means of] the 
implicit Euler method which leads to the non-linear equation 



where x*' is the solution at time t'' for k = 0, 1 , etc and h is the step t''*^ - 1^ This non- 
linear equation can be solved, for example, iteratively by [means of] the Newton method 



for n = 0, 1 , etc. 

[0025] In the parallel calculation, the master then only calculates the unknowns "m", the 
unknowns "s " representing fixed areas for the master. To be able to calculate a solution of this, 
the slaves or slave partitions calculate the corrected values sf""^ in each iteration step of the 
master and report these to the master. For this purpose, the slaves must solve the system of 
non-linear equations. 



where fj and D| are the corresponding equations and matrices for the slave "i". This non- 
linear equation can also be calculated by [means of] the Newton method and not all iteration 
steps have to be necessarily performed to convergence. 

[0026] The results of the slaves are then inserted into the system equations of the master 
whereupon the master can calculate the expression m^^^ . 



f( 
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[0027] In a next step 7. a new value "Cnew" for "C" is determined by the master and the slaves 
can also make suggestions. To achieve the aim, namely a value for "C" tending to zero or, 
respectively, less than a predetermined value e, and thus a calculation of the original circuit, the 
choice of the new value "Cnew" for "C" is made dependent on the difficulty of calculating the 
preceding step with the value "Cow" for "C". 

[0028] In this context, various procedures are conceivable, for example the choice of "Cnew" in 
dependence on an analysis of the number of iteration steps which the master needed for 
calculating the non-linear system of equations for the preceding value "Cow" for "C": 

C / 
0^%, IS n < n. 



Cold isrij <n<n2 
C . • 2 is n, > n 

old ^ 



where "n" is the number of iterations of the master and "ni" "x\2 are the parameters 
predetermined by the user. Naturally, instead of halving or doubling "Coid", respectively, other 
strategies for reducing or enlarging "Coid" are also possible. Additionally, the number of iteration 
steps needed by the slaves or slave partitions for solving their non-linear system of equations 
can also be taken into consideration, for example through the choice of 



n = maxj^n, X ^ > ' slavel,. . X ' slavep j 



where ni. siavej is the number of iteration steps of slave j, during the i-th iteration of the 
master. It must be noted in this context that the choice of a value for "C" only affects the 
efficiency of the method not the operating point itself. 

[0029] After a value for "C" of less than or equal to a predetermined value e has been reached, 
the calculation is terminated at branch 5 and the value for C can be set to "0" in a last step. The 
operating point of the original circuit is thus found as result 8 and can be output via output units 
such as for example, a screen, printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 
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[0030] Due to the method [according to the inventionl described above , an operating point of a 
very large electrical circuit can be advantageously calculated in parallel by a multiplicity of 
computers or processors and the disadvantage of known parallel types of calculation, namely 
lack of convergence with unfavorable starting values can be avoided. Due to the parallel 
calculation, the expenditure for calculation with regard to iteration steps and setting a value for 
"C" step by step, can be kept within tolerable limits due to the distribution to a number of 
processors. 

[0031] The operating point thus calculated is then used as basis for further analyses, for 
example the alternating-current analysis of a circuit. 

[0032] The invention has been described in detail with particular reference to preferred 
Q embodiments thereof and examples, but it will be understood that variations and modifications 
□ can be effected within the spirit and scope of the invention. 
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